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Abstract 

While designed primarily for X-ray imaging applications, the MedipixS ASIC can 
also be used for charged-particle tracking. In this work, results from a beam test at 
the CERN SPS with irradiated and non-irradiated sensors are presented and shown 
to be in agreement with simulation, demonstrating the suitability of the MedipixS 
ASIC as a tool for characterising pixel sensors. 
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1 Introduction 


As part of the upgrade of the LHCb experiment, scheduled for the second long shutdown 
of the LHC in 2018/19, the present microstrip-based Vertex Locator (VELO) is foreseen 
to be replaced by a silicon hybrid pixel detector with an ASIC dubbed “VeloPix” which 
will be derived from the Medipix family of ASICs. Like its predecessors Medipix2, Timepix, 
MedipixS, and TimepixS, the VeloPix chip will feature a matrix of 256 x 256 square pixels 
with a pitch of 55 qm. Prior to the arrival of the TimepixS ASIC, Timepix and MedipixS 
were the most suitable devices for the qualihcation of prototype sensors for the VELO 
upgrade. Until the end of lifetime of the npgraded experiment, the pixels closest to the 
beam line (r = 5.1 mm) accumulate a fluence of up to 8 x 10^^ 1 MeV Ugq cm“^. Qualifying 
silicon sensors in terms of radiation hardness is therefore a key element of the VELO 
upgrade R&D programme. 

Timepix silicon detectors have been characterised extensively in terms of charged- 
particle tracking performance [^. Since the ASIC has per-pixel information on the 
collected charge in terms of a time-over-threshold (ToT) valne, a direct measnrement of 
the charge deposition spectrum is possible. 

MedipixS on the other hand is a pure counting chip which was designed primarily 
for photon imaging applications. In order to measure the deposited charge in the sensor 
one therefore has to resort to indirect methods (Section]^. Among the large-scale ASICs 
in the high-energy physics commnnity, MedipixS has been the hrst to be based on IBM 
ISO nm CMOS technology. Compared to the Timepix (which was fabricated in 250 nm 
technology), it is expected to be more radiation tolerant and thus lends itself for testing 
irradiated sensors. 

In this work, we report on measurements with irradiated and non-irradiated MedipixS 
assemblies carried out in 2012 at the H8 beamline of the CERN North Area facility, using 
positively charged hadrons with a momentnm of 180 GeV/c. These measnrements are 
intended to provide a validation of the chip functionality and performance complementary 
to characterisation measurements using photon sources. In addition, they also represent a 
hrst step towards a comprehensive evalnation of the radiation hardness of silicon pixel 
sensors with the “Medipix footprint” of 55 x 55 qm^ pixels. 

2 Setup 

The Timepix telescope, described in Ref. j^, was nsed for reconstrncting the tracks of 
particles crossing the MedipixS device under test (DUT). In order to minimise the pointing 
error, the DUT was placed in the centre of the telescope. Eor reading ont the MedipixS 
chip, the “Merlin” data acqnisition system developed at the Diamond Light Source 
facility was used. The synchronisation of Timepix and MedipixS is straightforward, as 
both ASICs work in a “camera-style” frame-based readont mode. An external circnit 
implemented using NIM modules was used for sending shutter opening and closing signals 
to the telescope and the device under test. The coincident bring of two scintillators located 
npstream and downstream of the telescope was nsed for connting the nnmber of beam 
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Figure 1: DUT mount and cooling setup. A Peltier cooler (not visible in the photograph) is 
clamped between the cooling block and the TPG sheet onto which the chip is glued. 


particles traversing the telescope and the dnration of a shntter was adjnsted snch that 
50 scintillator triggers were accnmnlated in one frame. Dnring one spill (9.6s) typically 
abont 500 frames were recorded. 

The temperatnre of the irradiated samples was controlled nsing a combination of 
thermoelectric and CO 2 cooling, the latter being provided by a portable cooling plant j^. 
As can be seen on the photograph in hgure the chip was glued onto a tape of Thermal 
Pyrolytic Graphite (TPG) the other end of which was attached to the cold side of a Peltier 
cooler. The hot side of the Peltier cooler was put in contact with an aluminium cooling 
block through which GO 2 was circulating. With the chip switched off, a temperature 
of approximately —20° G was reached, rising to —15° G with the chip in operation. The 
setup was placed in a light-tight aluminium case which was flushed with nitrogen. The 
non-irradiated assemblies were measured at room temperature. 

In all measurements discussed below, the DUTs were n-on-p silicon sensors bump- 
bonded to MedipixS.l ASIGs. The ASIGs were operated in single-pixel high-gain mode 
and only one threshold (DAG THL) was used. Before taking data, a threshold equalisation 
was performed using the front-end noise as reference. The equalisation procedure consists 
of optimising the DAG values of two global current sources (ThresholdN and DAGpixel) 
and the threshold adjustment bits of each pixel, such that the THL value corresponding to 
the noise floor lies within a certain window for all pixels. For the non-irradiated assemblies 
a target window of 5 < THL < 25 was used. 

Prior to the beam test, calibration measurements using testpulses were performed to 
determine the relation between THL DAG and injected charge. To verify the viability of 
the testpulse method, data were taken with a ^"^^Am source for the two non-irradiated 
assemblies discussed below, and - after applying the calibrations obtained from the testpulse 
scan - the signal peaks were found to match within 2% between the two assemblies. These 
measurements were however made with a different readout system and equalisation 
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mask than used in the testbeam (at that time the testpulse functionality was not yet 
implemented in Merlin). 

Where possibl^ calibration measurements using test pulses were later (after the beam 
test) made also using the Merlin system. These calibration curves are used in the following 
for a relative comparison of the signals measured with different sensors. 

3 Measurements 

For each track reconstructed in the telescope, the intercept with the DUT plane is calculated. 
In order to suppress fake tracks a requirement on the track quality is applied and the 
tracks are required to include hits on all telescope planes. If an unused cluster with a 
centre of gravity within a radius = 110 pm around the track intercept is found, it is 
associated to the track and tagged as used. In case of multiple candidate clusters, the 
closest one is selected. The hit efficiency (or cluster finding efficiency) e is then given by 
the fraction of tracks with an associated cluster on the DUT. The pointing resolution of 
the telescope (~ 1.5 pm in the present configuration) allows one to probe the hit efficiency 
as function of the track intercept within a pixel cell. In the analysis we divide the pixel 
cell in 9 X 9 bins and calculate separate efficiencies for each bin. 

The most probable value (MPV) of the deposited charge can be estimated by scan¬ 
ning the hit efficiency as function of threshold. Assuming that the distribution of the 
collected charge can be described by a Landau distribution fi convoluted with a Gaussian 
distribution fc, the hit efficiency as a function of the threshold Q is given by 

OO 

e{Q) = j dxfi G) fc (a;) ■ (1) 

Q 

By htting the measured efficiency with Eq. ([^, the MPV and width of the Landau 
distribution and the a of the Gaussian can be determined (the mean of the Gaussian is 
fixed to zero). This method requires that the entire charge deposited by a track is collected 
by a single pixel. We therefore use the hit efficiency in the central bin of the 9x9 matrix 
for this measurement. 

To determine the spatial resolution, the distributions of the residuals between the 
X, y coordinates of the track intercepts and the associated clusters are calculated. The 
standard deviation of the residual distribution is used as a resolution metric (the pointing 
error of the telescope represents only a small correction when subtracted in quadrature). 
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Figure 2: Hit efficiency (for tracks crossing the centre of a pixel cell) in non-irradiated MedipixS 
assemblies with 100 p-m thick n-on-p sensors as function of threshold, for (left) assembly W20_B6 
and (right) assembly W20_J9. The insets show the distribution of the charge deposition in units 
of THL. The error bars represent the statistical uncertainty of the measurements. For several 
points the error bars are smaller than the symbols. 


4 Results 

4.1 Non-irradiated assemblies 

The measurements before irradiation were carried out using n-on-p active-edge sensors 
with a nominal thickness of 100 pm manufactured by VTT|^ The MPV of the charge 
deposition spectrum is thus expected to be around 7000 — 7500 electrons. 

Two assemblies, W20_B6 and W20_J9, were tested with beam. After the beam test, 
part of the backside metallisation was removed from the sensor on W20_B6, and by 
injecting laser pulses the depletion voltage Vhep was determined to be approximately —15 V. 
In the beam test, both sensors were operated at a bias voltage of —60 V, and were oriented 
perpendicularly to the beam. For each point in the threshold scan, a data set comprising 
typically 1 — 2 x 10® reconstructed telescope tracks was recorded. 

From the hit efficiency in the centre of the pixel cell as function of threshold (figure]^, 
the most probable value of the charge deposition spectrum is determined to correspond to 
THL ~ 111.2 for assembly W20_B6 and THL ~ 71.6 for assembly W20_J9. The statistical 
errors of the fit values are given in figure]^ Uncertainties due to tracking cuts, alignment, 
non-linearity of the THL DAC give rise to a systematic error on the MPV of ~ ±2.5%. 
The MPVs in terms of THL DAC values differ signihcantly between the two devices. This 

^One of the assemblies (W20_B6) was accidentally damaged after the beam test. 

^VTT Technical Research Centre of Finland, Espoo, Finland 
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Figure 3: Hit efficiency in non-irradiated MedipixS assemblies with 100 pm thick sensors as 
function of the track intercept within a pixel cell, for (left) assembly W20_B6 at THL/MPV~ 0.45 
and (right) assembly W20_J9 at THL/MPV~ 0.56. 


can be attributed to non-optimised settings of the DAC “Vcas” (for voltage cascode), 
which sets the reference voltage for some transistors in the pixel circuit and, if not tuned 
properly, can have an impact on the operating point of the circuit. 

As discussed in section testpulse calibration measurements using the Merlin readout 
system could only be made for assembly W20_J9. To facilitate the comparison between the 
devices, the results discussed below are therefore presented as function of the threshold-to- 
signal ratio, i. e. the applied THL DAC normalised to the respective MPV. 

Figure (left) shows the hit efficiency as function of the track intercept within a pixel 
cell for the lowest threshold-to-signal ratio (THL/MPV~ 0.45) covered by the threshold 
scan. At this threshold - which is high compared to a typical operational threshold of 
1000 electrons - the detector can be seen to be fully efficient (e > 0.99), except at the 
corners. With increasing threshold - as illustrated in hgure (right) - a drop in efficiency 
becomes noticeable at the borders of the pixel cell. This is a consequence of charge sharing 
due to diffusion, as can be seen from hgure which shows the average cluster size as 
function of the track intercept within the pixel cell. 

Figure 1^ (left) shows the average cluster size as function THL/MPV. At the lowest 
threshold, an average cluster size of 1.103 ± 0.001 is found, with single-pixel clusters 
constituting ~ 91.3% and two-pixel clusters ~ 7.7% of all associated clusters. The cluster 
size as function of threshold/signal follows the same shape for both sensors, exhibiting a 
minimum around THL/MPV~ 0.9 and a subsequent maximum around THL/MPV~ 1.6. 
With increasing threshold, an increasing fraction of the observed clusters is produced by 
primary particles which suffer collisions with large energy loss. These collisions give rise to 
energetic electrons which further ionise along their path and produce electron-hole pairs 
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Figure 4: Average cluster size as function of the track intercept within a pixel cell for assembly 
W20_B6 (100 pm thick sensor) at THL/MPV~ 0.45. Multi-pixel clusters are found predominantly 
at the edges and corners of the pixel cell. 


away from the trajectory of the primary particle. 

The residual distributions (in the x direction) at low threshold for one-pixel and 
two-pixel clusters are shown in hgure Averaged over all cluster sizes, the standard 
deviation of the residual distribution is found to be ~ 15.7 pm in both x and y. As 
expected from the dominance of one-pixel clnsters, this valne is close to the binary limit 
(55 pm/v^ 15.9 pm). Fignrej^ (right) shows the a of the residnal distribntion in x as 
fnnction of threshold/signal. The resolntion can be seen to deteriorate signihcantly after 
the threshold crosses the MPV. 

Figure also includes results obtained with a Timepix ASIC bump-bonded to a 100 pm 
thick n-on-p sensor from the same batch. The Timepix assembly was measured in the 
same beam test campaign and was operated at a threshold of 1000 electrons and a bias 
voltage of —GOV. The data were taken in ToT mode, but the values shown in hgure 
were calculated with the ToT values set to one to mimic the MedipixS behaviour. The 
results for the Timepix assembly are in agreement with the extrapolated resnlts from the 
MedipixS assemblies. 

To nnderstand better the observed shapes of resolntion and cluster size as function of 
threshold, a simple simnlation nsing the Gar£eld-|-+ toolkit 10 was nsed. The primary 
ionisation process is calculated using the Heed program 11 , which in addition to the 
energy loss by the traversing charged particle also simnlates the ionisation cascade from 
high-energy (“delta”) electrons prodnced in the interactions of the charged particle with 
the silicon medinm as well as the spatial distribntion of the resnlting electron-hole pairs. 
Each electron is snbseqnently transported throngh the sensor, based on the drift velocity 
and diffusion coefficient as function of the electric held. A one-dimensional approximation 
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Figure 5: Average cluster size (left) and standard deviation of the x-residual distribution (right) 
at perpendicular track incidence as function of threshold-to-signal ratio for assemblies W20_B6 
(empty circles) and W20_J9 (full circles). Results for a Timepix assembly with the same type of 
sensor are shown with full squares. The error bars represent the statistical uncertainty of the 
measurements. For most points the error bars are smaller than the symbols. 


for the electric field is usec0 

As can be seen from figure the features in the measured cluster size and resolution 
as function of threshold-to-signal ratio are reproduced by the simulation, provided that 
the spatial extent of the ionisation pattern is taken into account. This corroborates the 
conclusion that these features are due to “delta” electrons. 


4.2 Irradiated assemblies 


A set of MedipixS.l assemblies, with VTT n-on-p active-edge sensors from the same batch 
as the non-irradiated sensors discussed above, were irradiated at the Ljubljana TRIG A 
reactor 12 to a 1 MeV neutron equivalent fluence of 5 x 10^^ cm“^ and subsequently 
annealed for 80 minutes at 60° C. From lab measurements with a laser, the effective 
depletion voltage Vdep of these sensors after irradiation was measured to be between —85 
and —105 V. One of the assemblies, W20_H5, was characterised in the beam test. The 
sensor, which has a pixel-to-edge distance of 100 pm, was operated at a bias voltage of 
— lOOV as operation at higher bias was inhibited by the onset of electrical breakdown 
(hgurej^. 

In general, the irradiated assemblies exhibited a higher dark count rate compared 
to the non-irradiated ones, which was attributed to electrons from the /3-decay of ^®^Ta 


^The electric field is assumed to vary linearly between E = (V — Vdep) /d at the sensor backside and 
E = {V + Vdep) /d at the implants, where d is the sensor thickness. 
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Figure 6: Residual distribution for (left) one-pixel and (right) two-pixel clusters extending over 
two columns in a non-irradiated 100 pm thick sensor (assembly W20_B6, at THL/MPV~0.45). 



THL / MPV THL / MPV 


Figure 7: Simulated average cluster size (left) and standard deviation of the x-residual distribution 
(right) at perpendicular track incidence as function of threshold-to-signal ratio. Full (empty) 
symbols show the results with (without) the spatial extent of the ionisation pattern being included 
in the simulation. 


produced by neutron activation during the irradiation of the ASICQ This - presumably in 
combination with other radiation effects - resulted in problems during the equalisation 

^ The presence of in the irradiated assemblies was confirmed by gamma spectroscopy measurements 
performed by the CERN radioprotection group. 
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Figure 8: Leakage current as function of reverse bias voltage at T = —13° C after irradiation to 
0.5 X 10^^ 1 MeV neq cm“^ (100 (J.m thick n-on-p active-edge sensor with 100 |um pixel-to-edge 
distance). 


procedure for some assemblies, as well as a larger threshold dispersion. The optimised 
values of the DACs which control the global currents (ThresholdN and DACpixel) were 
larger (by a factor 2-3) compared to the non-irradiated assemblies. 

As can be seen from hgure the MPV of the charge deposition spectrum corresponds 
to a THL DAC of ~ 53.9. After correcting for the differences in gain and offset using 
testpulse calibration curves, this value is found to be approximately 8.5% lower than the 
MPV of the non-irradiated sensors. 

As in the non-irradiated case, the cluster size spectrum is dominated by single-pixel 
clusters. Cluster size and a of the x-residual distribution as function of threshold-to-signal 


ratio (hgure 10) follow closely the corresponding curves for the non-irradiated assemblies 
(hgure 1^. 


Figure pT] shows the hit efhciency as function of the track intercept within a pixel cell 
at the lowest measured threshold (THL = 40). In the centre of the pixel cell, an efhciency 
of 0.93 ± 0.01 is found, compared to 0.97 ± 0.01 for the non-irradiated assemblies at the 
same threshold-to-signal ratio (~ 0.74). 

A further beam test measurement was performed with a Medipix3.1 ASIC bump- 
bonded to a 200 pm thick n-on-p sensor manufactured by CN]V0 The sensor featured 
two guard rings and was diced at a distance of 400 pm from the border of the pixel 
matrix. The assembly (W20_D6) was also irradiated at Ljubljana, but to a higher huence, 
2.5 X 10^® 1 MeV Ueq cm“^. During the equalisation procedure, a large threshold dispersion 
was observed, such that the upper limit of the equalisation target window needed to be 
increased to THL = 40. In addition, a signihcant fraction (~ 10%) of the pixels needed 


^Centro Nacional de Microelectronica, Barcelona, Spain 
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THL 

Figure 9: Hit efficiency in MedipixS assembly W20_H5 with 100 pm thick n-on-p sensor as 
function of threshold, after irradiation to 0.5 x lO^^lMeVngq cm“^. The inset shows the 
distribution of the charge deposition in units of THL (Landau distribution convoluted with a 
Gaussian). 


to be masked. Tracks crossing a masked pixel were thus excluded from the efficiency 
measurements for this device. 

Because of time constraints, data were taken only at three THL values with this device, 
such that determining the MPV of the charge deposit spectrum was not possible. 

Figure [T^ shows the efficiency as function of the track intercept. The dependency of 
the efficiency on the applied threshold and bias voltage is shown in figure 
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Figure 10; Average cluster size (left) and standard deviation of the x-residual distribution (right) 
at perpendicular track incidence as function of threshold-to-signal ratio for irradiated assembly 
W20_H5. 


-i I I > |- 

0.3 0.55 0.65 0.7 0.63 0.49 0.29 0 . 13 * 



0.36 0.55 0.69 0.68 0.69 0.57 0.28 0.1 1 
J_I_I_I_I _ I _i_ 


jc [|j.m] 


Figure 11: Hit efficiency in MedipixS assembly W20_H5 with 100 |J.m thick n-on-p sensor as 
function of the track intercept within a pixel cell at THL = 40 (threshold/signal~ 0.74), after 
irradiation to 0.5 x 10^® 1 MeV ngq cm“^. The drop in efficiency at the edges of the pixel cell is 
more pronounced than in figure mainly because of the higher threshold-to-signal ratio. 
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Figure 12: Hit efficiency in Medipix3 assembly W20_D6 with 200 pm thick n-on-p sensor (diced 
at 400 (J.m from the last pixel) as function of the track intercept within a pixel cell at THL = 50 
and 300 V bias voltage, after irradiation to 2.5 x 10^® 1 MeVneq cm“^. The asymmetry of the 
efficiency profile is attributed to a small inclination of the assembly with respect to the beam. 




THL 

Figure 13: Overall hit efficiency in irradiated assembly W20_D6 with 200 pm thick n-on-p sensor 
as function of the THL value for different bias voltages (stars: — lOOV, triangles: —200 V, circles: 
—300 V). The assembly was irradiated to a fluence of 2.5 x 10^^ 1 MeVneq cm“^. The lines are 
drawn to guide the eye. 
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5 Conclusions 

First results on the performance of hadron-irradiated pixel detectors with the “footprint” 
of the upgraded LHCb VELO (55 x 55 pixels) have been obtained from a beam test 
at the CERN SPS, using MedipixS.l ASICs bump-bonded to thin n-on-p silicon sensors. 
As the MedipixS ASICs do not provide a direct measurement of the collected charge, 
threshold scans were used for characterising the sensor. 

Reference measurements with non-irradiated MedipixS. 1 assemblies were made in the 
same testbeam campaign. The threshold DAC value (THE) corresponding to the most 
probable value of the charge deposition spectrum was determined from the hit efficiency 
as function of threshold. The measured cluster size and resolution as function of the 
THL/MPV ratio are consistent between the assemblies and in good agreement with 
simulations. 

After irradiation with reactor neutrons to a fluence of 0.5 x 10^^ 1 MeV neq cm“^, a 
larger threshold dispersion and an increased power consumption of the ASIC were observed. 
A direct comparison between the performance of irradiated and non-irradiated sensors is 
complicated due to the fact that the non-irradiated sensors were operated at overdepletion 
while the irradiated sensor (due to early electrical breakdown) needed to be operated close 
to depletion. The difference in operating conditions may explain the small drop in efhciency 
and MPV measured with the irradiated sensor. The shapes of average cluster size and 
resolution versus threshold are not changed signihcantly compared to the non-irradiated 
sensors. 

The studies presented in this paper are being continued using the TimepixS ASIC, 
which has recently become available. 
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